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Fermi surface of superconducting LaFePO determined from quantum oscillations 
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We report extensive measurements of quantum oscillations in the normal state of the Fe-based 
superconductor LaFePO, (T c ~6K) using low temperature torque magnetometry and transport in 
high static magnetic fields (45 T). We find that the Fermi surface is in broad agreement with the 
band-structure calculations with the quasiparticle mass enhanced by a factor ~2. The quasi-two 
dimensional Fermi surface consist of nearly-nested electron and hole pockets, suggesting proximity 
to a spin/charge density wave instability. 

PACS numbers: 71.18.+y, 74.25.Jb, 74.70.-b 
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The recent discovery of superconductivity in ferro- 
oxypnictides [l|, l2j , has generated huge interest as as an- 
other possible route towards achieving high T c supercon- 
ductivity. LaFePO was among the first Fe-based super- 
conductor to be discovered and has a transition temper- 
ature of up to T c w 7 K pi . This compound is isostruc- 
tural with LaFeAsO, which is non-superconducting and 
has a spin-density wave (SDW) ground state [J] , but be- 
comes a relatively high-T c superconductor (T c rs 25 K) 
when electron doped [3(. By changing the rare-earth ion, 
T c reaches 55 K in SmFeAsOi_ x F x [Ref. i]. Theoret- 
ical models suggest that the pairing mechanism in the 
Fe-based superconductors may be mediated by magnetic 
fluctuations due to the proximity to a SDW 0, H, 0|- 
Knowing the fine details of the Fermi surface topology, 
its tendency towards instabilities as well as the strength 
of the coupling of the quasiparticlcs to excitations is im- 
portant for understanding the superconductivity. 

Quantum oscillations provide a bulk probe of the elec- 
tronic structure, giving detailed information about the 
Fermi surface (FS) topology and mass renormalization. 
To observe quantum oscillations samples must be ex- 
tremely clean and the upper critical field must be low 
enough for the normal state to be accessed; LaFePO is 
a material which fulfils both these requirements. The 
tetragonal layered structure of LaFePO is made of alter- 
nating highly conductive FeP layers and poorly conduct- 
ing LaO layers stacked along the c axis |l|, hence the 
Fermi surface is expected to be quasi-two dimensional 
[8| . Here we report extensive measurements of quantum 
oscillations in torque and transport data. We find that 
the Fermi surface of LaFePO is composed of quasi two- 
dimensional nearly-nested electron and hole pockets with 
moderate enhancement of the quasiparticle masses. 

Single crystals of LaFePO, with dimensions up to 
0.2 x 0.2 x 0.04 mm 3 , and residual resistance ratios 
p(300 K)//o(10 K) up to 85, were grown from a tin flux 
9]. Single crystal x-ray diffraction gives lattice param- 
eters a = 3.941(2) A, c = 8.507(5) A, and La/P posi- 



tions z La = 0.148901(19), z P = 0.63477(10) in agreement 
with previous results [l| . Torque measurements were per- 
formed with piezoresistive microcantilevers [lfj] down to 
0.3 K on different single crystals from the same batch 
(T c w 6 K); one in Bristol up to 18 T (sample B) and 
another crystal at the NHMFL, Tallahassee, up to 45 T 
(sample A) . Interplanc electrical transport has been mea- 
sured on a third sample (sample C). 

Figure la shows raw torque measurements versus mag- 
netic field at low temperatures (T = 0.3 K) at various an- 
gles between the magnetic field direction and the c axis 
|12| . At low fields we observe behavior typical of a bulk 
anisotropic type-II superconductor in the vortex state 
[ll| . The signal is reversible, indicating weak pinning 
of vortices. The upper critical field (Fig. la) is strongly 
anisotropic, varying between [AqH C 2 = 7.2 T when the 
magnetic field is parallel to the ab plane and ^oH C 2 — 
0.68 T when B\\c. For magnetic fields above ~ 9 T, we 
observe oscillations periodic in inverse field, which arise 
from the de Haas- van Alphen (dHvA) effect. The oscilla- 
tory signal is clearly visible by subtracting a monotonic 
background (a fifth-order polynomial) from the torque 
data (dHvA in Fig. lb) or transport data [Shubnikov-de 
Haas effect (SdH) in Fig.lc]. 

From the fast Fourier transform (FFT) spectra of the 
oscillatory data (Figs. Id and le) we can identify the 
dHvA frequencies F, which are related to the extremal 
cross-sectional areas Ak , of the FS orbits via the Onsager 
relation, F = {Ti/2ne)Ak- From the evolution of these 
frequencies as the magnetic field is rotated from being 
along the c axis (6 = 0°) towards the a axis (8 = 90°) 
[121 ]. we can construct a detailed three dimensional pic- 
ture of the shape and size of the Fermi surface. Four 
frequencies have significantly larger amplitudes than the 
others in both samples A and B as shown in Figs. Id and 
2a. We label a\ and «2 the two closely split frequencies 
at F w 1 kT (AF « 35 T) and the two higher frequen- 
cies, f3\ and 02. Besides these intense features (and their 
harmonics (l0|), in sample A (which was measured to 



M»H(T) 




Sample C 



8 = 43° 



1 2 

F(kT) 



40 45 

HoH(T) 




qU ^cy (T) 10C 



FIG. 1: (color online) a) Torque measurements on LaFePO 
measured at T — 0.35 K and in magnetic fields up to 18 T 
for different magnetic field directions. The arrows indicate the 
position of H C 2- b) Oscillatory part of the torque (dHvA) and 
c) resistance (SdH) in magnetic fields up to 45 T. d) Fourier 
transform spectra (field range 15-45 T) for two different mag- 
netic field directions, e) The angle dependence of the Fourier 
transform spectra for the field range 10-18 T (sample B). 



much higher fields) we see several other frequencies with 
smaller amplitudes; these are labeled 7, S, e and have 
frequencies in the range 0.7-1.7 kT. An extremely weak 
signal, T), was observed only in sample B (measured up to 
18 T) which we believe originates from a small misaligned 
crystallite [12J. The observed frequencies correspond to 
a fraction varying between 2.8% to 9% of the basal plane 
area of the Brillouin zone. They are significantly larger 
than those observed in the double-layer Fe-As compound, 
SrFe2As2 [la ], which is non-superconducting, and is likely 
to have a reconstructed Fermi surface at low tempera- 
tures due to a spin-density-wave ground state. 

Figs, le and 2a show the angular dependence of the 
main frequencies and their amplitudes. For a purely 
two-dimensional cylindrical Fermi surface the dHvA fre- 
quency varies like l/cos# and deviations from this indi- 
cate the degree of warping for a quasi-two-dimensional 
cylinder. As shown in Fig. 2b, the orbits fa and fa both 
originate from sections of FS which have significant warp- 
ing, but with opposite curvature (i.e., maximal and min- 
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FIG. 2: (color online) a) Angle dependence of all observed 
frequencies. Different symbols correspond to sample A (open 
circles), sample B (filled circles) and sample C (diamonds). 
Possible harmonics of the main frequencies are shown by tri- 
angles and the dashed lines indicate their calculated location. 
Solid lines are fits to a 1/cos 9 dependence, b) Angle de- 
pendence of F cos 6. Solid lines are calculations for a simple 
cosine warped cylinder, c) The temperature dependence of 
the Fourier amplitude for 6 — 32° (data are offset for clarity). 
Solid lines are fits to the Lifshitz-Kosevich formula 15]. The 
obtained values of the effective mass are listed in Table I. 



imal areas respectively) . The angle dependence of these 
orbits is well described by Yamaji's analysis of simple 
cosine warping of a two-dimensional cylinder [14], as is 
the large increase in the amplitude of the oscillations at 
~ 45°, where the frequencies cross (see Fig. le). These 
observations indicate that fa and fa arise from the same 
FS sheet. For the a orbits the warping of the FS sheet 
is very small, but the almost identical amplitudes, effec- 
tive masses (see below) and frequencies strongly suggests 
that these orbits arise from a single, weakly warped, FS 
sheet. 

The effective mass of the quasiparticles on the various 
orbits were determined by fitting the temperature de- 
pendent amplitude of the oscillations to the conventional 
Lifshitz-Kosevich formula [15|. as shown in Fig. 2c. The 
obtained masses range between 1.7 m e and 2.1m e and are 
listed in Table I (m e is the free electron mass). 

We now compare our experimental observations with 
predictions of the density functional theory calculations 
of the electronic structure of LaFePO. Our calculations 
were made using the WIEN2K code with the experimen- 
tal lattice parameters and atomic positions and including 
spin-orbit interactions [16|. The resulting band struc- 
ture is in good agreement with that reported by Lebegue 
8]. The density of states at the Fermi level arc derived 
mainly from Fe and P bands suggesting that the carri- 
ers flow mainly in the 2D FeP layer. The Fermi surface 



TABLE I: Effective masses (m*) and frequencies from dHvA 
data for samples A and B. Calculated band masses (mb) in 
LaFePO for shifted and unshifted bands. Orbits are labeled 
by band number and the location of their center. 



Branch 


F(kT) 


m 


/rn e 


Orbit 


rrib/m 


e 


(exp.) 




Sample A 


Sample B 


(calc.) 


unshifted shifted 


Oil 


0.985(7) 


1.9(2) 


1.8(1) 


1Z 


1.0 


1.9 


a> 


1.025(7) 


1.9(2) 


1.8(1) 


2r 


0.9 


0.7 


fa 


1.91(1) 


1.7(2) 


1.8(1) 


2Z 


1.2 


1.1 


fc 


2.41(1) 


1.8(2) 


1.9(1) 


3r 


1.8 


1.1 


5 


1.36(2) 


1.8(3) 




3Z 


2.9 


2.5 


7 


0.73(2) 


1.7(3) 




4M 


0.7 


0.7 


c 


1.69(2) 


2.1(3) 




4A 
5M 
5A 


0.9 
0.8 
0.9 


0.9 
0.7 
0.8 



Not Shifted 






Shifted 


b 


band 


■ 


— 2 / 




— 4 / 




— 5 J 


v 




■ X. 


I , . i . , 



-45 45 90 -45 45 90 

6 (degrees) (degrees) 



c 



■A 
.M 



CD 

.r 

G 



J I 



«v 




\ 


I 


1 


-A 






1 \ 


■M 


) ( 


1 

-r 

1 




o © 




FIG. 3: (color online) a) Calculated dHvA frequencies ver- 
sus angle for unshifted bands compared to experimental data 
(shown in light blue) b) As in a) but for shifted bands (see 
text), c), d) Two-dimensional (110) cuts of the Fermi sur- 
face and e), f) two-dimensional cuts in the FXM plane (001) 
for unshifted and shifted bands, respectively. Solid black lines 
indicate the Brillouin zone. A possible inter-band nesting vec- 
tor, Q ~ [it, tv, 0], between electron and hole bands is shown. 



mainly comprises small, slightly warped tubular sections 
running along the c direction. There are two hole cylin- 
ders centered on the Brillouin zone centre (r) and two 
electron cylinders centered on the zone corner (M) (see 
Fig. 3c-f ) . In addition, there is a small three-dimensional 
(3D) hole pocket centered on Z. The spin-orbit interac- 
tion makes small but significant changes to the Fernri 
surfaces; most notably it breaks the degeneracy of the 
bands crossing the Fermi level along the XM line such 
that the two cylindrical Fermi surfaces no longer touch 
and it also increases the separation of the elliptical hole 
pocket and the two cylindrical hole surfaces (Fig. 3c). 

The frequencies of the extremal dHvA orbits obtained 
from the calculated band structure are compared with 
the experimental data in Fig. 3a. The calculation pre- 
dicts that there should be 9 frequencies (two for each tube 
plus one for the 3D hole pocket) in the range 0.7 — 3 kT 
(for 9 ~ 0°), which is broadly similar to what is observed 
experimentally. In particular orbits fa and fa closely re- 
semble those expected from the larger electron cylinder 
in both frequency and curvature. The shape and splitting 
of orbits u\ and a-i are similar to that of the smaller elec- 
tron cylinder. The larger amplitude of these oscillations 
indicates that they both suffer significantly less damping 
than the other orbits (we estimate a mean free path of 
« 1300 A and 800 A for a and (3 orbits, respectively). 
This may be a natural consequence of the fact that both 
of these electron orbits originate from the same piece of 
Fermi surface in the larger unfolded Brillouin zone (cor- 
responding to the Fe-sublattice |5j). Assuming the above 
assignment for the a and (3 frequencies, the three remain- 
ing frequencies (7, S, e) would then correspond to hole 
orbits, although their exact assignment is less clear. For 
these hole orbits the scattering is a factor ~ 2 larger than 
the electron orbits (Figs. Id and le). 

Performing small rigid shifts of the energies of the elec- 
tron and hole bands improves the agreement with exper- 
iment. The two bands giving rise to the electron orbits 
a and (3 are shifted by —83 meV and —30 meV, respec- 
tively and the hole bands by +53 meV (Fig. 3b). The 
band that gives rise to the 3D pocket (which we do not 
observe experimentally) influences significantly the de- 
gree of warping along the c-axis of the hole cylinders; if 
this band was absent we would have a better agreenrent 
between our data and calculations (Figs. 3a and 3b). 

A two-dimensional cut in the TXM plane for our cal- 
culated unshifted and shifted Fermi surface (see Figs. 3e 
and 3f) shows that the electron pockets at the corner of 
the Brillouin zone (M) have almost similar shapes and 
sizes to the hole pockets at the centre of the zone (r). 
Nesting requires a perfect nratch between the size and 
the Fermi surface topology of the electron and hole pock- 
ets and this could stabilize a spin density wave (magnetic 
order) or charge density wave (structural order) with a 
wave vector Q « [tt, tt, 0] @, 0, H- LaFePO is non- 
magnetic [l8| but importantly our results suggest that 
size and the shape of the electron and hole pockets are in- 
deed very similar (Figs. 3e and f) implying that LaFePO 



may be close to nesting and hence to a spin/charge den- 
sity wave instability. 

For an undoped LaFePO the volume of the electron 
and hole sheets should be equal (compensated metal) but 
shifting the bands (as described above) creates an imbal- 
ance of ~ 0.03 electrons per formula unit. A possible 
explanation for this imbalance can be related to a small 
amount of electron doping in our crystals due to oxy- 
gen non-stoichiometry (~ 1.5% oxygen deficiency which 
is below the resolution of our x-ray data) . 

The band structure calculation allows us to estimate 
the many-body (electron-phonon and electron-electron) 
enhancements of the quasiparticle masses over their band 
values. Table I shows the band masses for all the pre- 
dicted orbits as calculated, and after application of band 
shifts. For the identified bands a moderate renormal- 
isation is found, with m* /rnb = (1 + A) w 2. For 
the hole orbits the calculated frequencies do not corre- 
spond exactly to the experimental ones, however, as the 
band masses for most of the hole orbits are in the range 
0.7 — 1.1 m e and the measured m* values are rs 2m e , 
a similar level of enhancement is likely. Reported mea- 
surements of the electronic specific heat coefficient (7) 
vary between 7-12 mJ/mol K 2 [1, [II GJ]. In 2D ma- 
terials 7 can be estimated from the renormalised dHvA 
masses, summed over all of the observed FS sheets. Four 
quasi-2D FS sheets with an effective mass ~ 2m e would 
give 7dHvA ~ 6 mJ/mol K 2 , close to the lower end of the 
experimentally observed values, without including any 
contribution from a 3D pocket. 

A recent ARPES study of LaFePO [HI showed an 
electron-like FS pocket centered at the M point and two 
hole sheets at the T point. The electron sheet and the 
smaller hole orbit have areas close to the experimentally 
observed dHvA frequencies. When compared to band 
structure calculations, the bands found in the ARPES 
measurements are renormalized by a factor two, similar 
to the mass enhancements we have found. This suggests 
that the renormalization occurs over the whole band, 



rather than being localized to energies close to the Fermi 
level, as found in Sr2Ru04, [2Q|. One of the hole pock- 
ets seen in ARPES is much larger (> 12 kT) than any 
dHvA frequencies observed here; its presence creates a 
significant charge imbalance of one electron per unit cell, 
suggesting that this feature is related to surface effects 

m 

As in the other ferrooxypnictides, in LaFePO the band- 
structure is very sensitive to the position of the P atom 
[2ll |22j |. Considering the P atom in the position op- 
timized (the 'relaxed' structure) results in hole orbits 
which are much too small to explain our data. Much 
larger band energy shifts are needed to bring it into ap- 
proximate agreement and we find that the band structure 
calculated with the experimental P position provides a 
better description of the data. 

In conclusion, we have found experimentally that the 
Fermi surface of the superconductor LaFePO, is in broad 
agreement with band structure calculations with a mass 
enhancement of about 2 due to many-body interactions. 
The difference in amplitudes of the dHvA signal between 
the hole and electron pockets is an indication of different 
scattering rates affecting these orbits. The near-perfect 
matching between the hole and the electron orbits that 
we observe, suggests that LaFePO may be close to a 
spin/charge density wave transition and that magnetic 
fluctuations are an important ingredient in the physics 
of the Fe-based superconductors. 
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